Modeling flux noise in SQUIDs due to hyperfine interactions 
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Recent experiments implicate spins on the surface of metals as the source of flux noise in SQUIDs, 
and indicate that these spins are able to relax without conserving total magnetization. We present 
a model of 1// flux noise in which electron spins on the surface of metals can relax via hyperfine 
interactions. Our results indicate that flux noise would be significantly reduced in superconducting 
materials where the most abundant isotopes do not have nuclear moments such as zinc and lead. 
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Although there have been significant advances in su- 
perconducting qubits, they continue to be plagued by 
noise and decoherence. Low frequency 1// flux noise 
[H in superconducting quantum interference devices 
(SQUIDs) is one of the dominant sources of noise in su- 
perconducting flux 0, IH and phase [ij qubits. Recent 
experiments indicate this flux noise arises from the fluc- 
tuations of spins residingon the surface of normal metals 
[H and superconductors Q . These spins have a high den- 
sity (~ 5 x 10 17 m -2 ), and may arise from local electron 
moments in localized states at the metal-insulator inter- 
face 0- 

One early model of flux noise due to spins proposed 
that the spin of an electron in a surface trap is fixed, but 
that the orientation of the spin can change when the elec- 
tron hops to a different trap @. However, the density of 
defect traps needed to explain the experiments was or- 
ders of magnitude larger than what is estimated to exist 
in a typical glassy material [3] • Another model suggested 
that spin flips of paramagnetic dangling bonds occurred 
as a result of interactions with tunneling two-level sys- 
tems mediated by phonons [l]. However, to obtain 1/ / 
flux noise, the maximum two level system energy splitting 
would have to be a few mK which is orders of magnitude 
smaller than accepted values. 

There is some experimental indication of interactions 
between the spins Q leading to the theoretical sug- 
gestion that flux noise is the result of spin diffusion 
via Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tions between the spins 11 1. RKKY interactions 
between randomly placed spins produce spin glasses, and 
Monte Carlo simulations of Ising spin glass systems show 
that interacting spins produce 1/f flux and inductance 
noise 

E3- 

In addition, RKKY interactions conserve the total spin 
and magnetization, so the total magnetization should not 
change. However, the Stanford group measured the to- 
tal magnetization of small isolated gold rings and found 
that the total magnetization is not conserved since the 
magnetization follows the externally applied ac magnetic 
field [fj . While this does not rule out magnetization con- 
serving interactions such as RKKY, it does imply that 
the spins must (also) be involved in interactions that do 



not conserve total magnetization. (The angular momenta 
contained in the external magnetic field and in the elec- 
tric current induced in the gold ring are orders of magni- 
tude too small to conserve total angular momentum by 
accomodating the change in angular momentum of the 
spins associated with the change in total magnetization.) 

The 3 possible interactions that do not conserve to- 
tal magnetization are spin-orbit, magnetic dipole-dipolc, 
and hyperfine interactions. The orbital angular momen- 
tum of a neutral gold atom is zero, so we can ignore spin- 
orbit interactions. The dipole-dipole interaction between 
2 electrons is of order 1 mK if we use a separation of 1.4 
nm corresponding to a spin density of 5 x 10 17 vnT 2 . This 
is much smaller than the hyperfine contact interaction 
which is of order 70 mK in hydrogen, for example. This 
implies that the hyperfine interaction dominates. In sup- 
port of this is the fact that scanning SQUID microscope 
experiments [B[ found that the magnetic susceptibility 
of spins on silicon is 5 to 20 times smaller than that of 
metals and insulating A10 x . This is consistent with hy- 
perfine interactions since the only isotope of silicon that 
has a nuclear spin and hence hyperfine interactions is 
29 Si which has a natural isotopic abundance of 5%. Note 
that spin angular momentum is conserved in hyperfine 
interactions where there is a spin flip exchange between 
the nuclear spin and the electron spin, but the magnetic 
moment is not conserved since the gyromagnetic ratios 
of the nuclear and electron spins differ by about 3 orders 
of magnitude. Previous authors [S, 11| have pointed out 
that flux noise cannot be directly due to fluctuating nu- 
clear spins because the frequency range and magnitude 
of nuclear flux noise would be much lower than what is 
seen experimentally. However, this does not rule out the 
possibility that the electron spins that are responsible 
for flux noise can relax via hyperfine interactions with 
nearby nuclear spins. 

In this paper we present a model of flux noise in which 
electrons residing in harmonic traps undergo spin ex- 
change with nearby nuclei via the hyperfine contact inter- 
action. The relaxation time T\ of a given electron spin is 
dominated by exchange with the nearest nonzero nuclear 
moment. In materials where not all the isotopes have 
a nuclear moment, the distance to the nearest nucleus 
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with a nuclear spin could be quite large. For example, 
the only isotope of Pb with a nuclear moment is 209 Pb 
which has a 22% natural isotopic abundance. The dis- 
tribution of distances between trapped electrons and the 
nearest nucleus with a magnetic moment gives rise to a 
distribution of electron spin relaxation times T\, which 
in turn results in 1/f noise up to 10 MHz. Unlike the 
model [lH of spin diffusion via RKKY that found white 
noise at low frequencies in contradiction to experiment, 
we find that 1/f flux noise extends down to 10 -5 Hz. 

The Hamiltonian of an electron spin S that is in an 
external field H oxt and that has a contact hypcrfine cou- 
pling to nearby nuclear spins I, is given by [131 ] 

ft = Hq + T-Lhyp 

Ho = — .9M-B H cxt • S 
= -hw S z 



n), (i) 



where I z = ±1 and S z = ±1. We ignore the dipolar 
interaction between the electron and the nuclear spins 
because it is much smaller than the contact hypcrfine 
interaction 14J. We choose H ext parallel to z. fiQ is the 
permeability constant. The external field could be due 
to an applied external magnetic field or to the magnetic 
field produced by local electric currents, and r are the 
coordinates of the i-th nuclei and the electron, go is the 
free-electron g-factor and j n is the nuclear gyromagnetic 
ratio, /is is the Bohr magneton of the electron. Taking 
the expectation value of Whyp with respect to the electron 
wavefunction ip(r) yields 



n 



(2) 



where "0( r O is the wavefunction of the electron at the 
position of the ith nucleus. () e indicates the expectation 
value of the Hamiltonian "Hh yp with respect to the elec- 
tron wavefunction. This can be written in the standard 
form of a hypcrfine interaction: 



(H hyp ) e = £(4 f l ; ) • s 



(3) 



is the hyperfine coupling constant between the ith 
nuclear spin and the electron spin. A typical hyperfine 
frequency, e.g., for hydrogen, is fhj — A^f/h ~ 1.4 GHz 
where h is Planck's constant. We can also express this in 
terms of an effective random field produced by the 
ith nuclear spin: 



where 



(Whyp)e — /. .90MbHr ■ S 



Hk = TrM07«?#(ri)| 2 Ii 



(4) 



(5) 



To find the effective field on the electron due to the 
ith nucleus, we use Eq. ([5]) and assume that the electron 
is in a harmonic trap with a ground state wavefunction 



(6) 



where £ 2 = h/m e uj, m e is electron mass and fl is the 
frequency of the harmonic oscillator. (If we assume 
that the localized wavefunction decays exponentially as 
exp(— r/£i) where is the localization length, then we 
still obtain 1/f noise up to logarithmic corrections in the 
frequency.) This gives 
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(7) 



where Aq is a constant. 

The nuclear spin dynamics can be characterized by a 
correlation time tq that is roughly the time scale over 
which the nuclear spins keep their orientation 15 1 



q(t)^(H^(t + T)H^(t)) 

= (( J ff^ j ) 2 )cxp(-|r|/r ), 



(8) 



where j — x,y, z are the components of the random mag- 
netic fields. The Fourier transform of this correlation 
function is 
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1 + (27T/) 2 T 2 ' 



(9) 



We will regard tq as a constant. If tq is determined by 
spin diffusion via nuclear dipole-dipole interactions, then 
To ~ l/Dq 2 where D is the spin diffusion constant and 
q = tx I a where a is the typical distance between nuclear 
spins. 

The electron spin dynamics is given by the Bloch equa- 
tion [15l. 1 161] : 



d_ 
~dt 



Hb H cxt x(S) 
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(S x 



T, 
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(S y )y, 



(10) 

where T\ is the spin-lattice relaxation time and Ti is the 
spin-spin relaxation time defined by 15 , 1(| 



= ({(HU 2 ) + ((Hky) 2 )) J 
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(11) 

(12) 
(13) 



Here we used Eq. ©. 

Each nuclear spin produces an electron relaxation time 
Ti(ri) that we can obtain by plugging Eq. ([7]) into 
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Eq. (dJ) to yield: 



C 2 r 



Tx( ri ) e 4 (l+^ 2 ) 



(14) 



where C and 6 are constants. To find an effective relax- 
ation rate 1/Tf , we sum or integrate over all the nuclei 
that lie within the electron wavefunction. The relaxation 
rate will be dominated by the nuclear spin that is clos- 
est to the electron. Let r n denote the distance between 
the center of the electron wavefunction and nuclear spin 
closest to it. Then we obtain 



1 

Tf 



Y — 

4 r 1 



ri(r) 



7r6 

2^2' 



(15) 



where a is the lattice constant. We now need to average 
over all the localized electron spins on the surface. Since 
the electrons are uniformly and randomly distributed on 
the surface, r n has a distribution P(r n ). For a 2D square 
lattice, the distribution of r n is P(r n ) = 2Trr n /a 2 for 
< r n < a 1 2. The surface on which the spins sit can 
be disordered, possibly resulting in an exponent for r n 
that is different from unity. However, it is reasonable to 
assume P(r n ) — Airf^ 1 /a 1 where A\ is a constant and 
7 G (2, 4) (7 is of order the dimension). Since r n depends 
logarithmically on Tf for both an electron wavefunction 
in a harmonic trap and an exponentially localized elec- 
tron, the actual form of P(r n ) is not important. The 
resulting distribution of Tf is 



spin fluctuations. From the fluctuation-dissipation theo- 
rem, the low frequency (fko <C kT) spin noise is propor- 
tional to the imaginary part of the spin susceptibility that 
we can derive from the Bloch equations in Eq. (|T0|) fljj ]. 
The frequency dependence of the noise at low frequencies 
is determined by the z-component of the susceptibility 
17j . The resulting spin noise power is 18 1 

1/Tf 



<S z (w) 



dTf P{Tf)sech 



Arf GCh \k B T U ' 



k^T 



(1/Tf) 2 
(18) 



where the limits of integration have a wide range with 

u 1l,mm < 1 < w ?T,max- 

To relate S z (uj) to the flux noise, we need to know how 
a spin couples magnetically to the SQUID. The effective 
flux <I> c ff produced by the spin magnetization on a loop 
with current I is 11 



cff 



S(r)B(v) 



dr, 



(19) 



where S(r) is the surface spin density operator and B(r) 
denotes the probing magnetic field due to the current 
and, if applicable, an externally applied field. Consider a 
SQUID made from a strip conductor (where the width of 
the strip is d) circular in shape with radius R (measured 
from the center of the loop to the middle of the annulus) . 
(For a square SQUID with circumference L and width W, 
we replace R and d by £/27r and W, respectively.) If the 
penetration depth A is much smaller than the width, the 
current density at x near the center of the strip is J(x) — 
2J/M)[1 - {2x/ d) 2 ]- 1 ' 2 for (-d/2) + A < x < (d/2) - A 
[111 Il9f . This current density produces the magnetic field 
B(x) — fioJ(x)/2. Using this in Eq (fl"9|) . we obtain the 
flux autocorrelation function [111 120(| : 



P(Tf ) 



At 



27/2-1 



In 
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1 

Tf' 



(16) 



We can simplify the above formula by approximating the 
slowly varying function In (irbTf / 2£ 2 a 2 ) by its average 
value. Then the distribution function is inversely propor- 
tional to Tf, and we can write P{Tf) = D x /Tf where 
D\ is a normalization factor determined by 



dTfP{Tf) = 1 



(17) 



( T eS \ 
V l.max/ 



and (Timin) correspond to the minimum 
and maximum frequencies of the flux noise, and are de- 
termined by (r n ) max and (r„) m i n , respectively. Thus we 



find A - e 2 /2 [(r 2 ) max - (r 2 ) min ] = e/^r 2 

According to the Wiener-Khintchine theorem, the 
spectral density S(uj) of the noise is given by twice the 
Fourier transform of the autocorrelation function of the 



(*(t)*(Q)) 



i 2 



drdr' ( 5(r, t)B{r)S(r', 0)B(r'] 
(20) 



If we assume the spins are isolated, 

(S(r, t)S(r', 0)) = Q(^A - |r - r'\){S(r, t)S(r, 0)) 



= a 2 S z (t)Q(VA~\v-r'\) 

where a — l/A is the spin surface density, A is the av- 
erage area per spin, Q(x) is a step- function, and S z (t) is 
the spin fluctuation autocorrelation function. After inte- 
grating over r' (using J 9(vC4— |r — r'\)f(r')dr' « Af(r) 
for an arbitrary function /(r)), we obtain 



(*(t)*(0))=afe^|ln(| c j5,(t) 
The associated flux noise spectrum is 



= cr(g^sMo) 2 sech 2 



TlLO. 



e 



(21) 



1 



-1 (— 

k B Tj d U \2Xj Arl 2tt/' 
(22) 
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which give rise to 1// flux noise. Notice that the flux 
noise is proportional to the density a of electron spins. 
Note also that for materials with a low concentration of 
nuclei with magnetic moments, Ar 2 will be larger and the 
flux noise, which goes as 1/Ar 2 , will be smaller. This is 
consistent with the small susceptibility found on silicon 
samples Q where 29 Si is the only isotope with a nuclear 
moment and its natural isotopic abundance is only 5%. 

Let us estimate the flux noise magnitude at 1 Hz 
for a Josephson junction. We can set the tempera- 
ture factor to unity since (/itB-ffcxt = f^o) "C kgT 
for i/ext in the range of 1 to 100 G and T between 
25 mK and 10 K. Since (r„) m i n w and we esti- 
mate that (Vn)max/£ 2 > 30, we make the approximation 



that £ 2 /Ar 2 = 
R/d = 10, a 



e/2 [{rl\ 



5 x 10 



17,. 



,r 2 ) mi „] » 1/30. Using 
ln(d/2A) ~ 8.5, and 



g^ofJ-B ~ 11.3(/z$o)(nm), we estimate the amplitude of 



the flux noise to be S, 



1/2 
*,hf 



5 /xi'o/Hz 1 / 2 . This agrees 
with experimental values which are typically in the range 
of 1 to 10 /xfco/Hz 1 / 2 [l|2l|. 

Eq. (|22p gives the flux noise due to spins that are 
only on the surface of SQUIDs. However, paramag- 
netic spins have also been found on a diectric surface 
5]. So if unpaired spins also reside on the substrate, 
these fluctuating spins will also contribute to the flux 
noise, reducing the dependence on d. Let L be the self- 
inductance of the SQUID. Then we can follow Wellstood 
[22j and use the expression for the electromagnetic en- 
ergy E = LI 2 /2 = J \B(r)\ 2 d 2 r/(2fi ) to evaluate the 
integral in Eq. (J2DJ): 



<S*(/) = a^oig/J-B^sech 2 



Hujq 



k B Tj Ar 2 27rf 



(23) 



The dependence of the flux noise on the geometry and 
the substrate are included in L. This result agrees with 
recent experiments [23[ that found a nearly linear rela- 
tionship between flux noise and L when the inductance 
of the SQUID was enhanced by inductor coils. Further- 
more, our result also implies that flux noise and induc- 
tance noise should be correlated [24j. 

To summarize, we have presented a model of 1// flux 
noise in which electron spins on the surface of metals 
relax via hyperfine interactions. Since the electron spin 
relaxation time depends exponentially on the distance 
between the electron and the nuclear spin, the nearest 
nuclear spin dominates the spin relaxation process. The 
distribution of distances results in a distribution P(Tf ff ). 
Averaging over this distribution results in 1// flux noise. 
Experimentally, the SQUIDs producing flux noise are in 
steady state equilibrium, so the noise is normally station- 
ary and Gaussian. This is what we have assumed in our 
calculations. (Stationary means that the system, and 
hence the autocorrelation functions, are translationally 
invariant in time. For Gaussian processes, higher order 
correlation functions can be expressed as products of the 



two-point (lowest order) correlation functions 25[.) Since 
the magnetization sums over individual spins, the mag- 
netization noise, and hence the flux noise, is Gaussian if 
there are enough spins for the central limit theorem to 
apply. In both our model and experiment , non-Gaussian 
noise could arise in very small samples |26j . 

Our results indicate that flux noise would be signifi- 
cantly reduced in superconducting materials where the 
most abundant isotopes do not have nuclear moments 
such as zinc and lead. The only isotopes of zinc and lead 
that have nuclear moments are 67 Zn and 207 Pb which 
have natural isotopic abundances of 4% and 22%, re- 
spectively. Thus, compared to Nb SQUIDs, we would 
expect flux noise to be lower by roughly a factor of 25 
and 5 in Zn and Pb SQUIDs since the relevant factor is 
£ 2 /Ar 2 in Eq. ([22]). This is assuming that Nb, Zn, and 
Pb have approximately the same atomic arrangement on 
their surface with approximately the same density of sur- 
face spins. For experimentally relevant values, the flux 
noise expression in Eq. (|22|) does not have any tempera- 
ture dependence (sech(fiwo/kBT) s» 1 since fiwo "C fcsT). 
This is consistent with a quantum process such as hyper- 
fine exchange coupling, and with the plateau seen below 
0.5 K in plots of the flux noise versus temperature 21 1. 
The unusual temperature dependence of the flux noise 
that is experimentally found above 0.5 K 2l| may in- 
volve thermal fluctuations of the spins. 
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